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Abstract. The purpose of this paper is to derive rigorously the so called 
viscous shallow water equations given for instance page 958-959 in [A. Oron, 
S.H. Davis, S.G. Bankoff, Rev. Mod. Phys, 69 (1997), 931-980]. Such 
a system of equations is similar to compressible Navier-Stokes equations for 
a barotropic fluid with a non-constant viscosity. To do that, we consider a 
layer of incompressible and Newtonian fluid which is relatively thin, assuming 
no surface tension at the free surface. The motion of the fluid is described 
by 3d Navier-Stokes equations with constant viscosity and free surface. We 
prove that for a set of suitable initial data (asymptotically close to "shallow 
water initial data"), the Cauchy problem for these equations is well-posed, 
and the solution converges to the solution of viscous shallow water equations. 
More precisely, we build the solution of the full problem as a perturbation of 
the strong solution to the viscous shallow water equations. The method of 
proof is based on a Lagrangian change of variable that fixes the fluid domain 
and we have to prove the well-posedness in thin domains: we have to pay a 
special attention to constants in classical Sobolev inequalities and regularity 
in Stokes problem. 
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1 Introduction 

This paper deals with a free boundary problem of nonstationary Navier- 
Stokes equations in thin domain justifying the derivation of the viscous shal- 
low water equations without surface tension. For information, such a system 
may be used to describe long-scale evolution of thin liquid films : See Equa- 
tions (4.12a)-(4.12d) written in [TS]. Depending on the geometry, there are 
two kinds of free boundary problems for the motion of a viscous fluid. The 
first one describes the motion of an isolated mass of fluid bounded by a free 
boundary and the second one describes the motion of a fluid occupying a 
semi-infinite domain in R n bounded above by a free surface and below by a 
fixed part of the boundary. We will consider in the present paper the second 
problem and refer the interested reader to [TS], [J_2] for applications for thin 
liquid films and rivers. First mathematical results for free boundary prob- 
lems were local existence theorems: A first local existence theorem published 
in 1977 by V.A. Solonnikov, in [19] ; a local existence theorem for the 
second problem by J.T. Beale, in 1980, in [3J. The following years brought 
other local existence theorems for equations of motion of incompressible flu- 
ids: see for instance [JJ, [2], [22], [25], [26]. The next step in investigating free 
boundary problems for incompressible Navier-Stokes equations was to obtain 
global existence theorems initiated by J.T. Beale in his paper [4] published 
in 1984. This paper was devoted to the motion of a fluid contained in a three 
dimensional infinite ocean. The first global existence result concerning the 
motion of a fixed mass of a fluid bounded by a free surface appeared, in the 
paper [20], by V.A. Solonnikov in 1986. The method used to prove the 
global existence results in the two papers mentioned above are completely 
different. J.T. Beale examined the surface waves after transformating it 
to the equilibrium domain through an adequate change of variable (dilata- 
tion in the vertical variable) while V.A. Solonnikov applied Lagrangian 
coordinates and this way transformed the considered drop problem to the 
initial domain. Other global existence theorem can be found. All global 
existence results for incompressible fluids are obtained for initial data suf- 
ficiently close to an equilibrium solution. Note that A. Tani, in [22] [23J, 
used the strategy developped by V.A. Solonnikov for free surface prob- 
lem to get a local existence result through the use of Lagrangian coordinates 
to fix the domain and the method of successive approximations. In con- 
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trast to J.T. Beale, the regularity of the solution obtained by A. Tani 
is sharp. Note that in J.T. Beale existence of arbitrary T > but for 
sufficiently small initial data in dependence on T is also proved. The re- 
sult is obtained in the same class of functions than the local existence. To 
the author's knowledge, in all these studies Dirichlet boundary conditions 
at the bottom are considered instead of Navier-Slip boundary conditions. 
Only recent papers with some rigid-lid assumptions consider Navier-Stokes 
equations in thin domains with Navier-Slip conditions, see for instance [9], 
[14] . In this paper, we will adopt the Lagrangian transform strategy in order 
to prove existence and stability of the viscous shallow water solution. To 
do that, we have to prove the well-posedness in thin domains with a spe- 
cial attention to constants in classical Sobolev inequalities and regularity in 
Stokes problem. We also prove the stability result in the Lagrangian form 
as in the paper written by D. Hoff, see [13], where uniqueness of weak 
solutions of the Navier-Stokes equations of multidimensional, compressible 
flow is proved using the Lagrangian formulation. We remark also that, using 
the Lagrangian formulation, no surface tension term is needed compared to 
[8]. Remark that there exists, at least, two different shallow- water type sys- 
tems mainly derived from two different bottom boundary conditions choices 
for Navier-Stokes equations with free surface: friction boundary conditions, 
no-slip boundary conditions. The system based on friction boundary condi- 
tions assumption is studied in the present paper and was mentionned in |18j . 
The system based on no-slip assumption has been formally derived in [S] and 
recently justified in [8] assuming non zero surface tension. Reader interested 
by shallow-water equations and related topics introduction is refered to the 
recent Handbook [6]. 

The paper will be divided in four Sections. The first one presents the 
Navier-Stokes equations with free surface and the viscous shallow- water ap- 
proximation. The second section is devoted to the linearized Navier-Stokes 
equations around the approximate solution built from a solution of the vis- 
cous shallow-water equations. We give anisotropic estimates and Korn's in- 
equality in thin domains and then prove the well posedness of the linearized 
Navier-Stokes equations. The main theorem and its corollary are proved in 
Section 4 that means well posedness and convergence of the full Navier-Stokes 
equations in Lagrangian form. On the last section, we draw some conclusions 
on this paper and present some perspectives to this work. 
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2 Navier-Stokes equations with free surface 
and viscous shallow-water approximation 

In this section, we write 3c? Navier-Stokes equations for a thin layer of an 
incompressible Newtonian fluid with a free surface. We then show how 
to obtain an approximate solution of these equations in the shallow water 
scaling: for that purpose, we will use smooth solution (see W. Wang and 
C.J. Xu [27]) of a viscous shallow water model fl6]) derived for instance in 
A. Oron, S.H. Davis and S.G. Bankoff [H] (see also J.-F. Gerbeau 
and B. Perthame [12J). We construct formally a second order approxi- 
mate solution (with respect to the aspect ratio e which measures the relative 
shallowness of the fluid layer) which will help to solve the Cauchy problem 
for Navier-Stokes equations with free surface and provides also convergence 
to viscous shallow-water equations. For this purpose, we formulate the equa- 
tions in Lagrangian coordinates. 

2.1 The equations and an existence theorem for the 
associated viscous shallow-water system 

In this paper, we consider the motion of a thin layer of Newtonian fluid with 
a constant density, set to 1. The differential operators V^div^ will concern 
the horizontal coordinates x = (xi,^). In the following, the fluid velocity 
u is denoted (uh,Uv) where % G M 2 (resp Uy G M.) is the horizontal (resp. 
vertical) fluid velocity. Navier-Stokes equations are then written, in their 
non dimensional form as 

d t u + u.Vu + -j Vp = — j| + — div(D(«)) , 
^ ^ Ke 

divw = 0, (1) 
for all t > and (x, z) G Qt the fluid domain 

tt t = {(x,z) / x EX, z G (0, h(x, £))} , X = T n ,M n , n = l,2. 



The nondimensional numbers Fq, Re are respectively the Froude and Reynolds 
numbers defined as 

F n 2 = ^— , Re 



U 2 „ LU 



gL v 

with U a characteristic velocity of the fluid and L a characteristic wave- 
length. Navier-Stokes equations are completed by boundary conditions. The 
impermeability condition at the free surface reads 

d t h + u H \z=h-V x h = uv\ z =h, (2) 
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which means that the fluid velocity is parallel to the free surface. We also 
assume a Navier slip condition at the bottom 



u v \ z=0 = 0, d z u H \ z =o = lu H \z=o, (3) 
and at the free surface, we assume the continuity of the fluid stress 

Here D(u) = Vw + Va T is the total deformation tensor of the fluid and can 
be written as 

n / x ( D x (u H ) d z u H + V x u v \ , , , > T 

D ^ U) = { (d z U H + V X U V ) T 2d z U V ) ' D ^ H) = V ^ + (V^) . 

The vector n is the normal to the free surface: 



Since the fluid pressure is defined up to a constant, we suppose that p atm = 0. 
Note that equation (j3J) can be equivalently written as 



p\ z=h _ (2d z u v - {V x h) T D x {u H ){V x h)) 



— , (5) 



F 2 Re(l-|V x ./i| 2 ) 

a . „ | ^ 2d z u v -{V x h) T D x (u H )V x h\ 

o z u H + VxU V \z=h = \ D x (uh) 1 - |V -h\ 2 / \ z=h ^ x 

In this paper, we will assume that the fluid domain is relatively thin and set 
h = eh , Fq = eF 2 with e« 1. As a result, the Froude number F = -J=u 
with H = eL, the characteristic fluid height, is the physical Froude number 
as the ratio between fluid velocity and wavespeed of perturbations at the 
free surface. Moreover, following the method of derivation by Gerbeau and 
Perthame [12] to obtain viscous shallow water equations, we suppose that 
the friction at the bottom is small: 7 = £7. Let us define an approximate 
solution of Navier-Stokes system f Trf2"f3"fl|) . For that purpose, we introduce 
(h ,uo) a smooth solution of the shallow water system mentionned by A. 
Oron, S.H. Davis, S.G. Bankoff [18], see page 958-959 with C -> 00: 



d t h + div x (h m ) = 0, 
'2F 2 Re 



h 2 1 

h (d t u + uq.V x u ) + V x — = r—&iv x {hoD x (u )) (6) 



+-^-^x(h div x (u )) - 
He Re 
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For the moment, we carry out a formal computation: we will make more 
precise assumptions on this solution in the sequel. We present a second order 
apporximate solutions, that means we search for an approximate solution of 
ffTOl3lT4j) in the form h a = eh and 

z^ 1 z^ z^ z^ 1 

U a ,H = U +U 1 Z+U 2 —, U a y = K) +t«l Z+W 2 — +W z — , p a = p +Pl Z +P^ ■ 

In the main Theorem we need higher order approximate solution that can be 
obtained in the same way, see [8] for construction of high order approximation 
of Navier Stokes systems with a free surface. Inserting this ansatz in the 
divergence free condition, one has 

w\ = — div x u , w 2 = — div x ui, w 3 = — div 2 .M 2 . 

Moreover, the boundary conditions at the bottom imposes 

Wq = 0, U\ = S^Uq. 

Next, using the second equation of (pQ) and boundary condition ([5]), one 
obtains the approximation of p as 

2sF 2 , 2s 3 ^?F 2 , 

p a = eh - z — div^Mo div x u . 

Re Re 

Next, we insert u a and p a into the second equation of (JSJ): identifying order 
0(e) terms, one proves that 

u 2 + W x wi = (D x (u ) + 2dw x (u )Id)—^ - 

Remark. For the approximate solution, we can calculate an ansatz up to any 
order on derivatives of (ho,u ) and polynomial dependence with respect to z. 
We refer to [S] for more details. We end up at order which gives a remaining 
term of order (z 3 ) for interior equations and (e 3 ) at free surface. 

A Navier-Stokes type system is satisfied by the high order shallow water ap- 
proximation (u a ,p a ). The adequate ansatz satisfies 

8 t u a + u a .Vu a + ^Vp a = + ^div(D(u a )) + 0(z 3 ), 

div u a = 0(z 3 ). (7) 

Moreover, the boundary conditions at the bottom are exactly satisfied and 
the boundary conditions at the free surface read 

d t h a + u a ^ H \z=h a y hK - u a ,v\z=h a = 0(e 3 ), 
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For the sake of completeness, we recall the existence and uniqueness result 
related to the viscous shallow water system ([6]) that can be obtained using 
the same procedure than in [27] 

Theorem 1 Let s > 0, u (0), h o (0)-l G H S+2 (R 2 ), \\h (t = 0)-l|| H .+a( R a ) < 
1. Then there exists a positive time T, a unique solution (uq, ho) of such 
that 

M , fco-le L°°(0, T; H 2+S (R 2 )), Vw G L 2 (0, T; # 2+s (M 2 )). 

Furthermore, there exists a constant c such that if \\ho(t = 0) — l||H*+ 2 (R a ) + 
\\uo(t = 0)||h s + 2 (k 2 ) < c then we can choose T = +oo with the same control 
on (h ,u ). 

Remark: the smallness assumption on ho(t — 0) — 1 can be replaced by 
ho(t, x) > a > to prove local well posedness provided that we work in Besov 
spaces: see the paper of R. Danchin [TT] for more details. However, we 
will need here this particular assumption to prove the well posedness of free 
surface Navier-Stokes equations in dimension 3 and convergence to shallow 
water equations. In dimension 2 for free surface Navier-Stokes equations, 
we just need an assumption similar to the one of [IT] , namely that there is 
no vacuum. Of course the results are only local in time: for arbitrary time 
interval, we need a smallness assumption on h (t = 0) — 1 and uo(t = 0). 

For global existence of weak solutions to a related viscous shallow water 
system, a result has been proved in [7] assuming an extra turbulent drag 
term. 

Remark: In this paper, we will choose s > 4. Indeed, a careful analysis of 
the remainder shows that the contained fourth order derivatives of (ho,Uo) 
have to be Lipschitz. 

2.2 Reduction of Navier-Stokes equations to a fixed 
domain 

In what follows, we write the Navier-Stokes equations in a fixed domain: this 
is done using a Lagrangian formulation of the equations. We first search the 
fluid velocity and pressure in the form 

eF 2 „ 

U = U a + U, p = p a + -TZ—P 
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Navier- Stokes equations (DP) then reads 



d t u + (u a + u).Vu + u.V u a + 
divw = 0. 



Vp _ 1 
Re ~ Re 



div(D(u))+/ a , 



in the fluid domain fl et = {(x, z)/x G X, < z < e/i(a;,t)}. The boundary 
conditions are then written as 



My 1 2 =o = 0, <9 2 Wtf 



2 = 



£7^# U=o- 



At the free surface, one has 



(D(u) - V ) z=ef n e = g a , 
d t (eh) + (u ajH + u H )\ z=£h .V x (eh) 



u a ,v + u v \ 

Next, we introduce the Lagrangian change of variable 



z=eh- 



(9) 



(10) 



= u atH (t, X, Z) + u H (t, X,Z), ^ = u a ,v(t, X, Z) + u v (t, X, Z) . 

For the sake of simplicity, we assume that the initial fluid domain is Xx (0, e): 
the rigorous justification will remain true for | < h (t = 0) > — 1| <C 1 and 7 
sufficienty large. Denote x G X, z G (0,e), the coordinates in this domain 
and A the Jacobian matrix of the Lagrangian change of variable 

/ OX 



A = 

The chain rules are given by 



0: 



A 1 



d z 



dX \ 



V Xo m d Zo u 



y xo w T 



d Zo w 



V x w T 



d z u 
d z w 



A. 



The fluid height is defined implicitely as eh(t,X(t,x ,e)) = Z(t,x ,e) so 
that the impermeability condition is satisfied. Moreover, one has 



eV x h(t,X(t, ,x ,e)) 



,dX, 
dx ' 



V X0 Z(t,x ,e). 



The Lagrangian velocity is defined as u = u(t, X) and u defined on the 
fixed thin domain X x (0, e). We further introduce the Lagrangian pressure 
p = p(t,X). In that setting, Navier-Stokes equations are written as 



A~ A 

d t u + {A- 1 u).Vu a + 



Re Re 
div(v4~V) =a a , 



.{dw(A T V)-VA:V)+f a1 (11) 

(12) 
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where f a = f a (t,X), u a = u a (t,X), (VA : V)i = Y,j, k d k A i,Pj*- The 
matrix V is the transformed deformation tensor: 

V = {Vu)A- l A- T + A- T {Vu) T A- T . 
The boundary conditions at the free surface zq = e reads 

V\ zo=£ -p\ Zo=£ Id)n = 0, n= f -(& TV ^Z(t,x ,e) \ ^ 



whereas the boundary conditions at the bottom z = are written as 

,dX 
8xq 



dX 

u v \z Q = 0, det(— )\d zo u H \ zo=0 = £7Mh| 20 =o- (14) 



The main result of the paper is 

Theorem 2 There exists C, K > so that if we assume 1 1 ITq 1 1 2 < C 2 e 3 ^ 2 , 
then there is a unique strong solution of (Illj) - (ll4p such that 

sup \\u\\ 2 2 + / \\u\\l + \\d t u\\i + \\p\\l + \\d t p\\l < Ke 3/2 . 
(o,r) Jo 

This gives the expected existence result of Navier-Stokes equations with free 
surface namely the system ([T]), (j2J), © and its convergence to the solution 
of the shallow-water system in a Lagrangian form namely a control of the 
solution of flBjl- fTTO]) with respect to e. 



Corollary 1 Let {Hq,Uq) be the solution of the viscous shallow-water system 
fl§D such that (h — l,u ) small enough in C(0, T; H 2+S (R 2 ) with s > 4. Let 
(/iq,Wq) such that \\(uq — w )|t=o||2 ^ Ce 3 ^ 2 and h £ = h and satisfying 
the compatibility conditions. Then System (JTJ, ©, © is well posed and 
u e — u a satisfy in its Lagrangian form the estimates given by Theorem 2. 

Note that there is no assumption on the time interval: the solution of Navier- 
Stokes system is defined on the same existence time as the solution of the 
viscous shallow- water system (EJ). 

The associated linear system. Next we introduce a linear problem that will be 
usefull to prove the well-posedness of the full Navier-Stokes system. Indeed, 
this will be done using a fixed point argument. First, let us introduce the 
following Lagrangian system of coordinates 

^- = u (X ), = -Z div(wo)(Xo), 
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with the initial conditions X (0, .) = x , Z o (0, .) = z . It is easily proved 
that X , Z satisfy 

dX i dX \ 

= 0, Z (t,x ,z ) = z ho(t,Xo(t,xo)), det^-^-J =h^ 1 ( l t,X (t,x )). 

In what follows, we will consider a special linear problem that is obtained 
if we substitute X = X and Z = Z in Navier-Stokes equations written 
in Lagrangian coordinates and drop order one differential operators. More 
precisely, this system is written as 

^ + = lt dW{AToVo) + /; (15) 

with V = V(u)Aq 1 Aq t + Aq t {Vu) t Aq t . The diver gence condition reads 

div(A x m) = a. (16) 

The boundary conditions at the bottom are 

uv\z =o = 0, ho 1 d Zo u H \ Zo =a = £lu H \ Zo =o + 9i, (17) 
whereas the boundary conditions at the free surface read 

[AlV -pId)\ Z0=£ ^ = g 2 . (18) 

3 Linearized Navier-Stokes around the ap- 
proximate solution of viscous shallow- water 
equations 

This section is devoted to study the linearized Navier-Stokes equations around 
the approximate solution of viscous shallow-water fll5|16fl7|fTH|) that was in- 
troduced in the previous section. We begin by a subsection dealing with 
Sobolev inequalities, Korn inequality and regularity result for Stokes prob- 
lem in thin domain. In the second subsection, we give well posedness result 
for the full linearized system. 

3.1 Anisotropic estimates in thin domains 

In this part, we recall Sobolev inequalities that arises in thin domains and 
pay a special attention to the constants arising in such inequalities. We also 
consider a linear Stokes problem in a thin domain that will be useful to obtain 
estimates on the Sobolev norms of (u,p). 



10 



3.1.1 Sobolev inequalities 

In what follows, we use that mean of u\ in the vertical direction is to 
write Poincare, Agmon and Ladyzhenskaia anisotropic inequalities in thin 
domains. Here the domain Q is Q = X™ x (0, e). In this paper, we will use 
the following notations: = f° r an y ^ G N with the convention 

H°(tt) = L 2 (n). Similarly, one denotes — I ^ (xv 

Proposition 1 There exists a constant C independent of e such that the 
following Ladyzhenskaia and Agmon estimates hold true: 

IMU 6 < C £_3 ||w||i, IMU^ < Ce~2 ||m|| 2 . 

The proof of this proposition is found in [21]. From this estimate and the 
inequality ||wv||l2 < \\u\\l6 \\v\\l6, one easily proves the following tame es- 
timates for the products of functions and composition of functions. 

Lemma 1 We will use the following estimates: if t > | 

\uv\hs(X) < C (\u\Ht(X)\v\H»(X) + M# s (X)Mtf*(X)), 

C 

\\uv\\ H s(n) < —^(WuWmi^WvWH"^) + ||w||fr«(n) 11^11^(0)), 



We have also the following property: For all F : MP — > M. smooth enough with F(0) 
then 

\\F(u)\\ s <C(\\u\\ L ~)\\u\\ s . 
We will also need the following trace theorems: 



Lemma 2 Let u G H 1 is so that u\ z= q = then one has \u\i < C||it||i, with 

2 



C 

C independent of e. In the general case, one has \u\i < — ^= J|^ || i 



This result is easily proved with Fourier series, that reduces the problem to 
a Id problem. The trace of a product is estimated using the following result. 

Lemma 3 Let u G and v G H 2 then, one has 

C 

kit; I i < —= \u _i \\v k- 
v^ 1 2" 



2 



The proof of this lemma is found in [10] in a fixed domain: we combine this 
proof and proposition [1] to obtain the estimate of lemma [3j 
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3.1.2 Korn inequality in thin domain 

In what follows, we prove the following result. 

Proposition 2 There is a constant C > independent of e such that 

2 \\D(u)\\ 2 L 2 (n) + e^\u H \ z=0 \ 2 L 2 iX) > C\\u\\ 2 H i {n) , 
for all u G H l {Vi) so that divu = and My 1 2=0 — 0. 

Proof. This is a consequence of Poincare inequality and trace estimates 
that the L 2 norm of uh and Uy are estimated as 

lkvlU 2 (n) < Ce\\d z u v \\ L 2 {n) , \\u H \\ L 2 {n) < C(e\u H \ z=0 \ L 2 {x) +e\\d z u H \\ L 2 {Q) ). 

We will prove a refined inequality on the deformation tensor when X = T n 

nD{u)\\ 2 mn) >C\\Vu\\ 2 mn) . 

For that purpose, we decompose the functions into Fourier series in the hor- 
izontal variables, (uh,v>v) = ^2 ke i 2 ( u k{ z )^ w k{ z )) exp(ik.x). Note that a 
similar result is obtained when X = IR n , using Fourier transform. We find, 
using the divergence free condition 

2||£>(u)||!2 ( n) = 2^ / kl\u kA \ 2 + k 2 2 \u ky2 \ 2 + \k.u k \ 2 
kei 2 ^° 

re rz pz 

+ ^ / \k 2 u ki i + hu k:2 \ 2 + \u' kl + / hk.Ukl 2 + \u' k2 + k 2 / k.u k \, 
Jo Jo Jo 

\\^ u \\l 2 (q)=^2 / \k\ 2 \u k \ 2 + \k.u k \ 2 + \ / k.u k \ 2 + \u' k \ 2 . 

where u'(z) = Setting U k = f* u k , there exists C k G (0,2) so that 

f 2k\\U[\ 2 + 2k 2 2 \U' 2 \ 2 + 2\k.U'\ + \k 2 U[ + k l U' 2 \ 2 + \U'{ + k x k.XJ\ 2 
Jo 

+ f \U' 2 ' + k 2 k.U\ 2 > C k f \k\ 2 \U'\ 2 + \k.U'\ 2 + \k.U\ 2 + \U"\ 2 , (19) 
Jo Jo 

for all U G i/ 1 (0,e) 2 so that U\ z= q = 0. We prove the inequality for each 
modes. We must prove that there is C > independent of k so that C k > C, 
forall k G Z 2 . If = 0, the two quantities are equal and one can choose 
Co = 1. In what follows, we assume that k 7^ 0. Denote C k G (0,2) the 
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optimal value for the mode k. Using test function in T>(0, e), one can prove 
that, if Cfc ^ 1, an optimal solution necessarily satisfies the differential system 



D^ + kjDU! = k ± k 2 DU 2 , D 2 U 2 + k\DU 2 = k x k 2 DU x , G (0,e). (20) 

where D is the differential operator D = — \ k\ 2 Id. This space of solutions 
is 8-dimensional and reduced to a 6-dimensional one, using the condition 
U\ z=0 = 0. Moreover, the eigenfunction U can be written U(z) = U(\k\z), 
where U G <S, the vector space defined as 

U(z) — ^ {a]_z + bi sinh(z) + c^cosh^) — 1)) 

+ (a 2 sinh(2;) + b 2 z sinh(z) + c 2 z cosh(z)) . (21) 

where we have denoted k as k = \k\(c,s) and Oj,6j,Cj G R, % — 1,2. Then 
we have to show that there is C independent of k,e so that Q 2 (U,k) > 
C Qi(U, k), for all f/ G S and Qi the quadratic forms defined as 



I'M 

Q 2 (U, k)= 2c 2 (U[) 2 + 2s 2 {U' 2 ) 2 + 2{cU[ + sU' 2 ) 2 
Jo 

+ / {sU[ + cU' 2 f + {U'l + c(cC/i + sU 2 )) 2 + (U% + s(cC/i + sU 2 )) 2 , 
Jo 

I'M 

Q^U, k)= \Uf + (cU[ + sU' 2 ) 2 + \U"\ 2 + (cC/i + sf/ 2 ) 2 , 

with M = e\k\. Let us denote qi(k),i = 1,2 the symmetric matrix associated 
to the quadratic forms Qi(., k), % — 1, 2. If is the minimum value for which 
the inequality Q 2 (U, k) > Ck Qi(U, k) for all U is satisfied, then det(q 2 (k) — 
CkQi{k)) = and one is left with an eigenvalue problem. Let us study the 
roots of P(X) = det(q 2 (k) — X q\{k)). First, an integration by parts in the 
quadratic form Q 2 (., k) yields Q 2 (U, k) = Q 2 (U, k) + Qi(U, k) 

Q 2 (U, k) = 2(cU 1 + sU 2 ) (M) (cU[ + sU' 2 ) (M). 

The quadratic form Q 2 (.,k) is a product of two linear forms Q 2 (U,k) = 
(vi(k),a)(v2(k),a) with (v 1 (k),v 2 (k)) G M 6 and a G M 6 represents the coor- 
dinates of U in the basis of S. Then the matrix q 2 (k) associated to Q 2 (., k) is 
given by q 2 (k)a = (vi(k), a)v 2 {k) + (v 2 (k), a)vi(k). As a consequence, for any 
a G (vi(k), v 2 (k)) ± , one has q 2 (k)a = q\(k)a. One can prove that Vi(k),v 2 (k) 
are not colinear so that X = 1 is a root of P with multiplicity 4. There re- 
mains two extremal values which correspond to minimal and maximal value 
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of Cfc. We know that G (0,2): let us show that X = 2 is exactly the 
maximum value for C^. It is sufficient to show that there is a divergence 
free vector fields (u, w) so that J{u,w) = J(u,w) T and w(0) = 0, J(u,w) 
being the Jacobian matrix of u, w. This is done choosing a potential flow 
(u, w) = (Vt/>,<9 2 V0 so that i[) z (0) = and Atp = 0. Split the functions into 
their Fourier modes, one shows that it suffices to choose ip k {z) = a k cosh(|fc|2;) 
with a G / 2 (Z 2 ) chosen so that the H 1 regularity of u,w is satisfied. As a 
result, the polynom P is factorized in the form 

P(X) = T k (X - 1) 4 (X - 2)(X - A(k)) 

and there is a unique minimal value A(/c) which depends continuously of 
M = e\k\ and a = c + is G S 1 . We analyse A(k) when M ~ e and M — > oo. 
We have calculated with a formal computation software the expansion of P 
with respect to M at M = 0: 

P(X) = M 4 -{Z - lf(Z - 2) + C(M 6 ), 
9 

so that A(fc) ~ 1 as e — > and lies in a fixed compact set. Next we deal 
with the limit M — > oo: Assume that there is {7m so that MQi(Um, k) — 1 
and MQ 2 (Um, k) — > 0. This is equivalent to the existence of £7w so that 

1 M 2 (d7 liA f + ^ 2 , M ) 2 + + I^mI 2 + M,M + ^2,m) 2 = 

lim (c^, M (l) + sU 2 , M {l)) ° U[ ' M + M SU ' 2M = -1. 

Denote /m = cU\m + sU<i,m'- it is easily proved that /m(1) and f' M (l) are 
estimated as 

/m(1) 2 = 2^ f M (z)f M (z)dz<^, 

^/m(1) 2 < ^ j\f M {z)) 2 dz + ^[f z IfMf'Mldsdz < ± 

This yields the contradiction with liniM-s^oo ^~Vm(1)/m(1) = ~ 1 an d con- 
cludes the proof of the Korn's inequality. □ 

3.2 Well posedness of linearized Navier-Stokes equa- 
tions 

3.2.1 A special divergence free linear problem 

In this section, we consider the special linear problem obtained from f)15H18p . 
assuming h = 1, u = 0: the more general case \ho — 1| -C 1 and |« | <C 1 
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will be treated in the next section through a perturbation argument. In the 
case considered here, A = Id and the linear problem is 

V» 1 

d t u + = — div D{u) + f, divu = 0, (22) 
tie Re 

with the boundary conditions at the bottom 

uv\z =o = 0, d zo u H \ zo=0 = £7 u H \z =o + 9i, (23) 
and at the free surface 

(d Zo U H + ^x U V )\z =e = 92,1, (2d ZQ U V - p)\ Z0=£ = 92,2- (24) 

We prove the existence of a weak solution. The weak formulation of (I22|23|24p 
is written as 

d t J{% 4>) + J D(u) : D{4>) + g J(uhU=o, <Ph\z =o) = 

{f, <t>) + ^ J (92, 0U=e) - (.91, 0h|« o =o), (25) 

for all G (if x (fi)) 3 so that div0 = and 4>v\ Zo =o = 0. In what follows, we 
prove the existence (and uniqueness) of a solution of (1231) . 



Proposition 3 Let f G L 2 ((0, T), and G L 2 ((0, r),i?-5(T)). 

For any v,o G L^^fi), t/iere is a unique weak solution u G £ 2 ((0, T), ff 1 ^)) 
of (TJ3p and <9 4 m G L 2 ((0, T), (i7 1 (0)) / ) . Moreover, it satisfies the energy 
estimates 



max( 0)t ) IMIo + C / ||u||i< IKHo + 
Jo 

C I \\ff {Hmy + \92,2\\ + -(\92,i\\ + \9i\\). (26) 

JO 2 £ 2 2 

Proof. The proof follows from a classical Galerkin approximation process: 

we only prove here the energy estimate. Substituting = u in the weak 
formulation (1251). one obtains 



d\\u\\l 
dt 2 



- / D{u) : D{u) + e*f [ \u H \ 
* Jn Jt 



|2 

z=Olo 



{f,u) {HlyHl + {g 2 ,u\ z=e ) H _i H i -( gi ,u H \ z=0 ) H _i H i. (27) 
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Using the Korn's inequality on the functional space (H^ iv (Q)) so that 
uv\ z =o — 0, one proves that there exists a constant C > independent 
of £ so that 

Nl? < C( I D{u) : D{u) + £7 / \uh\ z =o 
v Jn Jt 

Moreover, using lemma [U one has the trace estimates 

\u v \z=e\i < C\\u v \\i, >/e\u H \i < C||w#||i. 

2 2 

Using a Young's inequality, one can prove that for any r] > 0, there exists 
c(rj) > so that 

\{92,l,U H \z=e) ~ {gi,U H \z=o) H -i H ^\ < ^-(\92,l\ 2 _i + | ^1 1 — i ) + 7?ll W #lll> 

l(«2^«v)| H -4 (fl .il < c{r])\g 2 , 2 \ 2 _ i + 

Inserting this inequality into (12 7p gives the estimate 

^\\u\\l + C\\u\\l < c(\\f\\l Hly + \g 2A \\ 2 + i(|, 2il | 2 + l^l 2 )). 

Then integrating this equation on the interval (0, t) yields (1261) . From the 
weak formulation, one easily proves 

r* i 

|^|| 2 ^iy < |K||2 + / + l#2,l| 2 1 + -(\92,l\ 2 _ 1 + \9l\ 2 _l) 

JO 2 £ 2 2 

Using an interpolation argument, one has u G C ((0, t), Lj iv ) and the unique- 
ness follows, letting / = g { = and u = into (l2"B"j) . This completes the 
proof of the proposition. □ 

This is a standard argument that there exists a pressure p e L 2 ((0, t), L 2 ) 

so that / p(.,t) = and 
Jn 

/ (u t , 4>) + D(u) : -D(0) - / pdiv0 + £7 / Wfr|z=o <M*=o 
Jn Jn Jx 

= f(fA)+ [ (92, <P\z=e) - (gi, <Ph\z=o) V0 G #\ 0^U=o = 0. (28) 
Jn Jx 

Moreover, the pressure term p satisfies 

/ |b||o < c(K||o+ / \\f\\ 2 {H iy + \92,2\-i + -{\92,i\ii + \9i\ 2 L 

JO X JO 2 £ 2 

where C is a constant independent of e. Next, we prove a result on the 
regularity of the weak solution. 
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Proposition 4 Let f G L 2 ((0, T), H 1 ) such that d t f G L 2 ((0, T), (tf 1 )') 

and gi G £ 2 ((0, T),H^), d t gi G £ 2 ((0, T),H~^) which satisfies the compat- 
ibility conditions 

(d zo u H + V xo u v ) zo=£ (t = 0, .) = 52,1(0, .), 
(<9 Z() W//U=o - £^u H \ ZQ=0 ){t = 0, .) = 51(0, .)• 

and divu = 0, Mo £ -^ 2 - T/ien i/ie weaA; solution u verifies 

max( 0) t)(||w||i + lkt||o) + ^ ||u(s, .) \\\ + K(s, .) \\\ds < 

ll«o||l + ||/(0, ||<fc(0, .)||! 

2 

/"* 1 

+ / 11/11? + ||/*||(Jfi)' + l<?2, 2 |! + \92,2,t\\ + -{\9iM + I^mI'i) 
JO 22522 

Proof. Let us first prove that there exists an "initial" pressure p G i? 1 (Q), 
satisfying the mixed Dirichlet/Neuman problem 

-Ap = -div(A Mo + /(0)), 

Po\z =e = 2d Z0 U y\ Z0=E - 52,2(0), 

<9 2o PoU,=o = £7div Xo M ,^| 2 =o + diva. ^i(0) + M0). 

This problem is well posed since /(0) G L 2 and ^(0) G i?^ with a com- 
pactness argument respectively on / G L 2 (0, t; H 1 ) fl -£^((0, t); (H 1 )') and 
g G L 2 ((0, £); n H 1 ^, t; H~%). We split this problem into three parts: 
first, it is easily proved using a Lax Milgram argument that there exists a 
unique solution p ,i £ -f^ 1 (^) so that 

-Ap ,i = -div(Aw + /(0)), po,i|*o=e = 9 z po,i|ao=o = 0. 

This is a straighforward computation (with Fourier series) to prove that there 
exists C, independent of e, such that 

lbo,i||i<c(||/(o)|| + |Kll2). 

Next we consider the problem 

Ap ,2 = 0, Po,2|*= £ = ^, ^p , 2 U=o = 0, 

with h = 2d z u y\ Zo=£ —52,2(0) G . Using Fourier series, this is equivalent 
to find Pq 2 so that 

d 2 d 
~d^ pk °' 2 + \ k \ 2pk °> 2 = °' Po, 2 l*= £ = h k , ^0,2^=0 = 0. 
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The solution is easily calculated and satisfies 
coshl&U , ,, d 



1 f y 

with C(y) = 5— / cos\i{2x)dx < Ci,Vy > so that llf 2 II 1 < Ci|^|i- 

cosh y J ' 2 

Similarly, the solution p ,3 of 

Ap ,3 = 0,Po,3\z =e = 0, <9 2 p ,3| 2 =o = e7div xo u 0l # | zo=o +div Xo 0i(O)+/v (0) | 2o=0 , 



satisfies ||po,3||i < Ci(||^o||i + |<?i(0)|i + ||/2(0)||o- Let us define an initial 
condition u G L 2 (f2) ( that plays the role of d t u(0)) as 

«o = -Vp + Am + /(0). 

This satisfies the following estimate 

IMIo + lbo||i < Ci(||«o|| 2 + ||/(0)|| + |^i(0)^ 2 (0)| 1), 

with C\ independent of e. We apply the previous proposition with this 
initial condition that lies in the space of initial conditions of weak solutions 
{u G L 2 / divu = 0, uv\ zo =o — 0} and f,g are replaced by (ft,9t)- There 
exists a unique weak solution u G L 2 (0, t; if 1 (f2)) associated to this problem 
and that satisfies the energy estimate 

max (0lt) ||«||o+ / Mi<C 1 ((\\u \\l + \\f(0)\\ 2 +\(g l ,g 2 )(0)\l 
Jo v 2 

ft 1 

+ / ll/t||(Hiy + b2,2 lt |ii + -(W 1 +\92,i,tti)- ( 29 ) 
Jo 2 £ 2 

Associated to this weak solution is the pressure pressure term p G L 2 (which 
plays the role of p t ) that satisfies the estimate L 2 (0, t; L 2 ) 



* PI! 2 <C 1 ((HH2 + ||/(0)||g + \9i(0),g 2 (^l 

ft 1 

+ / ||/t||(Hiy + l^| 2 _i + -(\9i,tU + l^, M | 2 _i). (30) 
Jo E 

Then, let us consider Ui — u + / u, p\ = p + / p. Integrating on 

Jo _ Jo 

the time intervals (0, t) the weak formulation on u and using the definition 

of po, one can prove that Ui is a weak solution of the initial problem with 
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ui(t = 0) = u . Then one finds that u — u±, p — p\ and u t = u, p t = p 
satisfies the estimate (T2"9~j) . This concludes the proof of the proposition. □ 

Finally, we check the space regularity of the weak solution. One can prove 
the proposition 

Proposition 5 Let f G L 2 (0, T; H 1 ) n F x (0, T; (H 1 )'), 

^GL 2 (0, T; Hl)nH\0, T; iH) 
and Uo G i? 2 satisfying the compatibility conditions 

d ZQ u H + V Xo My| 2o=e (t = 0) = g 2 ,i(0), (d Zo u H - £ju H )\z =o(t = 0) = ^(0). 
Then there exists a constant C independent of e such that the weak solution 
u G L 2 (0, T; iJ 3 (fi))niJ 1 (0, T; H\tt)),p G L 2 (0, T; J ff 2 (fi))ni/ 1 (0, T; L 2 (fi)) 
satisfies 



max( ,r)||«||I+ / \\u\\l + \\p\\ 2 2 + \\u t \\ 2 + \\p t \\ 
Jo 



<C 1 [\\u U + \\f(0)\\ 2 o+\(9i(0),92(0))\ 

r-T 



+ [ Wf\\l+\\ft\\ 2 {H lY + \g2,2\l + \92,2,t\-l 

Jo 2 2 

1 f T 

+- / (\gi,g2,i\l + \gi,t,g2,i,t\ti)- ( 31 ) 
£ Jo 2 2 



3.2.2 The full linear problem 

In this part, we consider the full problem with a divergence term and ho^l. 
Let us first deal with the case Hq — 1,Uq — (steady state) and consider the 
linear problem (I22|I23II24I) with the divergence free condition replaced by 

divu = h, (32) 

with h G L 2 (0, T; H 2 ) D i7 1 (0,T; L 2 ). In order to prove the well posedness 
of (I22ll23|l24p with the divergence condition f[3"2"j) . we will apply proposition HJ 
for that purpose, we introduce Uh = V0 so that 

A0 = h, —\ z= o = —\ z=£ = 0. 
oz oz 



19 



One can prove, with Fourier series, that there exists a constant C > inde- 
pendent of e such that 



+ Uth<c(\\h\\ 2 + \\h t \\ ). 

Next, we consider /, such that f = f + Au h — u^t and 

01 = 01 + £JUh\z=0, 01,2 = 01,2 - %d xz (j)\ z=e , 02,2 = 02,2 + 2<9 zz </>. 

This is a straightfoward computation to prove that there exists a constant C 
independent of e such that 



i<^(ll/l|i + IN| 2 + INIo), 

|01,01,2|| + |01,t,01,2,tUi < C (|01,01,2 1 § + |01,i,01,2,t| 2 I + ||/l||2 + K||o): 
|02,2|| + |02,2,t|-l < C(|02,2|| + l02,2,t|_i + -^(||/i|| 2 + K||o)J- 



The estimate on ||/i||(#i)' is obtained using the linear Navier-Stokes equa- 
tions: more precisely, one has ft — ft + ^u^t ~ u h,tt and 

dxvu h ,tt = h tt = divuu = div(/ t + — div(D(u t )) - 

Then, one obtains the estimate of Uh,tt (and thus f t ) in L 2 ((0, T), (H 1 )') norm, 
using elliptic regularity. Next, we apply proposition @] on the divergence free 
problem for the fluid velocity u — Uh and pressure p. Finally, we have proved 
the well-posedness of the linear problem f )22|23f24p with f |32|) and the solution 
(u,p) satisfies the energy estimates 



u \\l + IKIli + IMI2 + Ibtllo 



max(o,r)(||w||a + l|wt|lo) + / 

Jo 

<cf|KI| 2 + 11/(0)11* + K01, 02) (o)i! 



+ r 11/11? +ll/t||(Hiy + 102,211 + 102,2,^! 
Jo 2 2 

1 f T 

+ - / (INI2 + Il^llo+l0i>02,i|l + l0i,t,02,i,t|ii)- (33) 
£ Jo 2 2 

Let us now deal with the full problem (ll5|ll6|ll7|[T5j) . For the sake of simplicity, 
we assume a = 0, the inhomogeneous case a ^ is treated as previously. In 
order to deal with the more classical divergence free condition, we introduce 
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the Jacobian decomposition of the fluid velocity u = A u. The linear problem 
then reads (we also drop a linear term in the form L Q (u, Vw, x, t)) 

A T ,d t (A u) + ^ Vp = ^div(Po) + /, (34) 

with V = AIA V(u)Aq x Aq t + (Vm) t . The diver eence condition reads 

divw = 0. (35) 
The boundary conditions at the bottom are written as 

u v \z =o = 0, Hq 1 d zo u H \ ZQ=Q = e*/u H \ Zo=0 + ffi, (36) 
whereas the boundary conditions at the free surface read 

(v -pld^j | Z0=e e 3 = g 2 . (37) 

Similarly to the special linear problem treated in the previous section, one 
can obtain a weak formulation of ( I34||3"T|) : 

d t J (A u, A <p) +1-1^:^ + ^1 h \^u H \ Zo = , |^hL=o) 

= / (A 0tt (p,A u) + (f, 4>) + — / (g 2 , <f) H \z =e) ~ K X {^-9ii ^-^(Ph\z=o)- 
Jn Ke J x oxo oxq 

for all <p £ H 1 (Q) so that div0 = and 4>v\z =o- One prove the following 
theorem 

Theorem 3 There is rji > 0, r] 2 > such that for any (h ,uo) £ W 2, °° be a 
smooth solution of (0|) which first satisfies \ho(t — 0) — < r\\ and either 

• (small data, large time) Wt £ (0,T), \\h (t, .) — 1, u (t, .)||oo < or > 

• (large data, small time) Tsup tg ( 0T ) \\h (t, •) — l,Wo(t, . ) 1 1 v^ 2 -°° < ^2 
T/ien /or a// source terms f,gi, i = 1, 2 and /i so £/ia£ 

/ £ L 2 ((0, T); ^(O)) n ^((0, T); (tf 1 )'^)), 
£ L 2 ((0, T); #§(X)) n ^((0, T); iH(X)), 
h £ L 2 ((0, T); # 2 (fi)) n ^((0, T); L 2 (fi)). 

T/iere exzsfo a unique weak solution (u,p) of k22\23\2J$M) and this solution 
satisfies energy estimates [33) . 
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Proof. The proof of this theorem follows the one of proposition HJ one first 
constructs an approximate solution using a Galerkin method and then prove 
uniform estimates on this solution. The only difference here is that we have 
to prove the Korn's inequality when (h ,uo) ~ (1,0). Let us first assume 
that ho(t = 0) = 1. In the case of small data, one can write that 

/ W : V = I [ D{u) : D(u) + 0( m )\\Vu\\ 2 , 
Jn „ 1 Jn 

P s~i ~V P 

/ ^I-q-^hU^oI 2 = (l + ^Ofe)) / \u H \z =o\ 2 

Jx OXq J w 

Then, one can use the Korn's inequality proved in the case ho — 1 to show 
that for i]2 sufficiently small, there is C > independent of e so that 

Jn Jx ax o 

From the Korn inequality, we deduce energy estimates on the approximate 
solutions: this proves the existence of a weak solution. The proof of regularity 
is similar to the case of ho = 1, uq = 0: indeed the main task was to prove 
a Korn's inequality in this setting, which is done mainly using the fact that 
the Lagrangian change of variable is close to identity uniformly with respect 
to e. The case of large data and small time is simply treated by considering 
that the Lagrangian change of variable is close to identity for T small enough 
and one can prove a similar Korn's inequality in this case. The more general 
case h (t — 0) ^ 1 is treated by perturbation, simply through the change of 
variable (x,z) h-> (x,zh (t = 0,x)). □ 

4 Well posedness and convergence of the full 
Navier- Stokes equations in Lagrangian form. 

We write an iterative scheme to prove the well posedness of free surface Navier 
Stokes equations with Navier slip condition at the bottom together with the 
convergence of this set of equations to a viscous shallow water model. We 
will use Tychonov fixed point theorem, as in [TU], with now the functional 
space 

X = {(u,p) e L 2 ((0,T);if 3 x H 2 ) / \\(u,p)\\ x < C(\\u \\ 2 + £ 3 / 2 )}, 

II (Mil! = su p IMI2 + / \Hl + \\d t u\\i + \\p\\l + \\d tP \\ 2 . 

(0,T) Jo 
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The iterative scheme is written as 
d t u k ^ + A^u k +Wu a (X ) + A ° 7 

JrC e 

= ^(div(AK +1 ) - VA : V k+1 ) + F\ 

div (A^u k+1 ) = H k , 
with 

= a a (X , Z ) + (<7 a (X fe , Z k ) - a a (X , Z )) + dw((A 1 - A k l )u k ) 

and 

F k = (A l - Al 1 )u k .Vu a {X k , Z k ) + A^u k .V(u a (X , Z ) - u a (X k , Z k )) 
+ ^((V ~ + {Af " A o T )(div(^P fc ) - VA k : V k )) 

+ -^A - T (VA : V k - VA k : V k + div(A k T V k - A T V k )) 

+ (fa(X k , Z k ) - f a (X , Z )) + f a (X , Z ). 
This set of equation is completed with boundary conditions at the bottom: 



< +1 U=o = 0, ^< +1 Uo = e-fu h H +1 \ z=0 + G\, (38) 



with G\ = g 7 ^ ( detffix" ) — 1) u h)U=o an d boundary conditions at the free 
svirf ctcG 

(AK +1 -/ +1 /rf)| 2=£ e 3 = GL (39) 

where 

G\ = ^((^(X^Z fc )-^(Xo,Zo))+^(Xo,Zo) + (P fe -P fc )n+(P fe -p fc /rf)(no-n 
with 

V k = {Vu k )A k l Af + A- T (Vu k ) T A- T , 
V k = (Vu k )A^A T + Af(yu k fAf. 

The function u k+1 has the initial condition u k+1 \ t =o = u . Using the well 
posedness results obtained in the previous section, this set of equation has a 
unique solution and the weak continuity of the map is done. We first prove 
that (u k ,p k ) is bounded, uniformly on the time interval (0, T) and uniformly 
with respect to e, in a suitable norm provided that the initial condition is 
sufficiently small. The proof is done in 5 steps: estimates of the Lagrangian 
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coordinates (more precisely the deviation from the "shallow water" coor- 
dinates), the boundary terms, the divergence term and nonlinear terms in 
momentum equations. 



Step 1: Estimates of Lagrangian coordinates. Let us recall that Lagrangian 
coordinates satisfy the set of equations 

j(X k - X ) = u a , H (X k , Z k ) - u a , H (X , Z ) + u k H , 

-(Z k - Z ) = U a y(X k , Z k ) - U a y(X , Z ) + U k y. 

Let us first compute the basic L 2 estimate 

j t (\\x k - x \\l + \\z k - z \\l) < C(\\u a (t, .)\\ w ^)(\\x k - x Q \\l + \\z k - z \\ 

+ ||n fc || (||X fe -X || +||Z A; -Z || ). 
This is a straightforward consequence of Gronwall inequality that 

sup (\\X k -X \\l + \\Z k -Z \\l) <Te c ^ f T \\u k (t,.)\\ 2 dt, 
te(o,T) Jo 

with C = C(sup fg ( 0)T )(||w a ||v^i,°°)). There is no difficulty to extend this 
estimate to higher order space derivatives: 

sup (\\X k -X \\l + \\Z k -Z \\l)<Te c ^ f T \\u k {t,.)\\ldt, 
te(o,T) Jo 

with C = C(sup tg ( 0T )(||M a ||i4/i,oo)). Next, we compute boundary terms at 
the bottom: 

Step 2: Estimates of nonlinear terms at the bottom. In order to estimate the 
boundary term at the bottom, we will use lemma [2] and [21 The nonlinear 
term at the bottom reads 

G k = ^4U=o(^||||^ - l)U=o = eiu k H \ z= ,Hd x X k - d x X ). 
We first compute an estimate for |<7*|3. Let us fix to > 1 so as H to (T 2 ) C 
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L°°(T 2 ). We then have 



\Gl\§ <ej(\u k H \ z=0 \ to \F(d x X k - d x X )\s + \u k H \ z=0 \ 3 _\F(d x X k - d x X )\ to ) 
<eiC{\d x X k - d x XoU{\u k H \ z=0 \t \d x X k - d x X \ 



I to 

+\u k H \z = Q\s\d x x k -d x x \ t0 ) 
H\\ t0+ \\\x k -x,h , IKIl2ll xfc - x ol| to+ | 



u 



<e-fC(\X k - .Yoit 0+ i)( - "JT*"" - ' + 

'* + t \ flL.fcll II vk 



e y/e y/i yfi 

<^'<- ^ 1 )(ll" /S |L + i|l^ fc -^0||3 



||X fc -X | 



+ \\u k \\ 2 \\X k -X \\ to+i ). (40) 
Then it is easily seen, choosing to — §, that 

|G^(t,.)r3^<7^(Apil^l|L 2 ((o,T),H3))T sup (||« fe (t,.)ll2) / T ||n fc (t,.)||^. 
o 2 V £ *e(o,T) Jo 



Using lemma [3j we next compute an estimate of \d t g k \_i- One has 



d t G k = £7(^^| z=0 J-(9 :c X fc -^.X )+n H | 2=0 9 1 ^(^X fc -^X )(^(X fc -X ; 

Then, the following inequality holds: 

Ad t u k \_i 

\d t G k \_i<e^--^\\T(d x X k - d x X )\\ 2 

+\ UH \ z=0 g(d x x k - d x x )(d x u k H + Ji^x* - d x x )\_i 

<C7(y|||ri*|| La(( o in j I »))(||«*|| 1 + \\d t u%)\\X k -X \\ 3 

+^\d x u k \_i\\F(d x X k -d x X )u H \\ 2 

<C{^\\u k \\ L . mnm) ) ((\\u% + ||<Vl|i)||X fc - Xo\\ 3 + \\u%\\X k 

So that we obtain the estimate 

f \d t G k \\ < c(y|ii M fc n L2((0 , T) ^3 ) )r( r \\d t u k w\ + \\u k wi) 2 . 

Jo 2 v S J 

Next we compute interior estimates of "divergence" terms: this is done using 
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proposition [TJ and tame estimates of lemma [TJ 

Step 3: Estimate of divergence terms. One part of the divergence term 
H k = Tr((A k 1 - A 1 )\/u k ) satisfies 

\\H k || 2 < -^C(\\d x X k - d x X \\ oo )\\X k - Xoh\\ Vu%, 
As a consequence, one has 

't\\„m 



|fl*(f,.)|frft<— (/ \\u k \\l) 2 . 



Next, we compute an estimate of Ho^-fff ||o- On the one hand, one has 

WTidA, 1 - A^ l )Vd t u k ) || < ~^C(\d x X k - d x X U\\X k - X || 3 || Vila, 

V £ 

and 

\\Tr((d t A k -d t A )Vu k )\\ < llV^IUIIC^^-^^llo, 
\\Tr((d t A k - d t A )Vu k ) ||o < -Lt^ C (\d x X k - a^loc 

+ \d x u k \ oo )(\\X k -X \\ 1 + \\u k \\ 1 ). 

So that we obtain 

f \\WX<-C(J-\\A\»wn#>)) / PtuXdt \\u k (t,.)\\ldt. 



This concludes the estimate of nonlinear divergence terms. 

Step 4 : Estimates of "interior" terms. The interior terms are estimated 
with proposition [1] and lemma [TJ We first compute an estimate of ||(v4o T — 
A^ T )Vp k \\i. One has 

||(V - -4* Will < IK T - Af\U\p k h + ||V(V - A- k T )Vp% 

<^-AfU\P%^^ A7 ' T - Af)hn 



C{\d x X k -d x X 



< _^l^l \\x k - X \\ 3 \\p k \\ 2 . (42) 
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As a consequence, we obtain the estimate 



Jo £ \l £ 



u k \\ L 2 mn m)) I K(f,.)lls* / \\p%- 
o Jo 



Similarly, the time derivative of this term satisfies 



\\dt{(A T - A k T )Vp k )\\ 2 {H1 y 

<-C(^\\u k \\ L , mnH3) ) f \\u\t,.)\\tdt [ T \\d t p k (t,.)\\ldt 

e V £ Jo Jo 

The other nonlinear term with high derivatives that we consider is 

f k = Afdiv(A T k V k ) - A^ T dW(A^V k ). 
This is a lengthly but straightforward computation to show that 

T \\f%+w k \\** y 

< jC(^\\u k \\ L2mnH , } ) £ \\u k \\l(l Q T \\u k \\l+\\dtu k \ 

Step 5: estimates of boundary term at the free surface. We have to deal with 
the nonlinear term 

G\ = Al{{V k - V k )n + (V k - p k )(n - n)) \ z=£ . 
On the one hand, one has 

\(V k - V k )n\ z=E \s < C{^\\u k \\ L , mnH3) )\V k - V k \s. 
Moreover, this is a straightforward computation to show that 

\<D k <T> k \ nt H^ fc ~ l|A fc - Xo|| 3 ||M fc ||3 

\V -V \z<C{ -= ) . 
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As a consequence, one obtains 




For the term (V k — p k )(n — n Q )\ z=£ , we first remark that (n — no) 1 2=0 = 0. 
As a consequence, we obtain the estimate 

/ T |(P fc -/)(n-n )U|i 
Jo 2 




Hence, using the estimates on the linear problem and using that remaining 
terms in the Navier-Stokes type system (j7|) satisfied by the high order shallow 
water approximation (u a ,p a ) satisfy 

£(\\f a (X ,Z )\\j+\\dJ a (X ,Z )^ 
£ Jo v 

+ \gl t (X ,Z ),d t gl 1 (X ,Z )\ 2 _ 1/2 ^ < Ce\ 
We have then proved the following proposition. 

Proposition 6 Assume that \\u0W2 < > then for e small enough, one 

has, for any k > 1 

\\(u k ,p k )\\ 2 x <C(\\u \\l + e 3 ). 

At this step, there are two strategies: Either we use a Banach fixed point 
argument and in this case we need estimates for higher Sobolev norms that 
are obtained in a way similar to our previous calculation. In this case, we 
can prove that the sequence (u k ,p k ) is a Cauchy sequence respectively in the 
space C(0,T;H 2+r '(tt)) and L 2 ((0, T), H 2+1) {VL)) with < 77 < 1/2, see [3] for 
further references. In this case, the result is not sharp. To obtain the well 
posedness in critical spaces, one has to use Tychonoff fixed point procedure 
as in [TU]. There is no difficulty here to follow both strategies to obtain the 
well posedness result in our situation since the estimate in high norms is 
proved. 
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5 Conclusion 



In this paper, we have proved the well-posedness of free surface Navier-Stokes 
equations with a slip condition at the bottom and convergence to a viscous 
shallow water model in the shallow water scaling e — > 0, e being the aspect 
ratio. In contrast to our mathematical justification of a shallow water model 
for a fluid flowing down an inclined plane, we dropped the tension surface 
effect (which is important in order to deal with realistic situations). This is 
due to the fact that we worked in the more suitable Lagrangian coordinates. 
There is no difficulty to extend this result for fluids flowing down an inclined 
plane when the slope 9 is asymptotically small, 9 ~ e and when the uniform 
flow is stable (see [5] for more details). 

When the uniform flow is unstable, it is known that instabilities so called 
roll-waves appear: these are periodic travelling waves. Small amplitude roll- 
waves are proved to exist both in viscous shallow water equations [15] and 
free surface Navier-Stokes equations [T7] . However, the result of Nishida et 
al. is proved in the presence of surface tension whereas roll- waves are proved 
to exist in viscous shallow water equations without surface tension and are 
of large amplitude when the viscosity is small. In order to obtain a larger 
range of validity of shallow water equations, this would be of interest to prove 
the existence of roll-waves in Navier-Stokes equations without surface tension 
that are close to "shallow water roll-waves" in the shallow water scaling. 
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